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We investigate colierent propagation through a large optical density Tm''"*':YAG crystal. Using 
an ultra-stable laser, fiber filtering and site selection, we investigate the transmitted pulse temporal 
profile. The plane wave condition is satisfied by selection of the illuminated spot central area. 
We pay special attention to 7r-pulse transmission in the prospect of implementing optical quantum 
storage protocols. 
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I. INTRODUCTION 

Coherent response to optical excitation has been inves- 
tigated for a long time. Actually, mere linear absorption 
of an optical beam through a material slab represents the 
most elementary manifestation of this feature. Indeed 
the transmitted beam results from the coherent combina- 
tion of the incident beam and the instantaneous material 
response. It has also been realized early that coherent 
optical response is related to the collective quantum ex- 
citation of an atomic ensemble. As pointed out by R. 
Dicke [H, all the molecules are interacting with a com- 
mon radiation field and hence cannot be treated as inde- 
pendent. However, a semiclassical treatment is generally 
adequate to describe the coherent response pj, espe- 
cially if the number of atoms per A-cube is small enough 
so that superfluorescence cannot start up on the experi- 
ment timescale J2| • In the framework of the semiclassical 
description, McCall and Hahn were able to derive an area 
theorem that accounts for the propagation of an arbi- 
trary intensity pulse through an arbitrary optical density 
slab 3, giving rise to the notion of Self-Induced Trans- 
parency (SIT). To be coherent, the process must take 
place on a time scale much shorter than the atomic su- 
perposition state lifetime. Experimental demonstration 
in the nanosecond range was performed at once in both 
solids and atomic vapors at room temperature [1, 0]. 
The McCall and Hahn theorem only predicts the pulse 
area behavior. For a more detailed description, especially 
to investigate the pulse deformation, one has to solve the 
Maxwell-Bloch (MB) equation numerically. 

The connection with solitons have stimulated the in- 
terest in SIT. It was shown [4] that a 27r-area pulse 
with hyperbolic secant temporal profile can propagate 
without alteration through an absorbing medium. The 
breakdown of large area pulses into 27r-area components 
was also predicted [1] and observed Q. A couple of 
decades later, the discovery of Electromagnetically In- 
duced Transparency (EIT) in ensembles of three-level 
atoms renewed the interest in the coherent response 
of optically thick media. Actually a variety of induced 
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transparency processes have been investigated in three- 
level systems 0, J!]. The interplay of EIT with SIT was 
also considered [9|. In recent years, effort has been di- 
rected to the ultra-short pulse regime, where the usual 
assumptions, such as the rotating wave and slowly vary- 
ing envelope approximations may not be valid [lO|. 

Rare-earth ion doped crystals (REIC) represent an ex- 
cellent testbed for coherent optical response investiga- 
tion. Most of the limitations pointed out by Slusher in 
his review paper on self-induced transparency [l3| are re- 
laxed in these materials. The superposition state lifetime 
may reach hundred of microseconds and is not limited by 
a transit time across the light beams. No diffusive mo- 
tion can affect the spatial phase of the atomic states. The 
inhomogeneous broadening can be considered as infinite. 
Generally, the level degeneracy is completely lifted by 
the crystal field, so that the ions behave as true two-level 
atoms on the propagation timescale. Other features de- 
pend on the specific host matrix. In Yttrium Aluminium 
Garnet (YAG, Y3AI5O12 ) for instance, the D2 symmetry 
of substitution sites forbids the existence of permanent 
dipole moments. As a consequence, relaxation processes 
such as optical excitation induced instantaneous spectral 
diffusion (1ST) are strongly reduced. In this matrix, ap- 
propriate crystal cutting and adequate polarization ori- 
entation make the ions in different sites interact with the 
same dipole moment projection on the incoming field pd| . 

Despite of such attractive features, very few SIT co- 
here rit p ropagation data have been collected in REIC 
Il4| . Ref [13] reports on experimental work in 
Tm''+:YAG. However coherent excitation conditions are 
not clearly satisfied since the laser is not stabilized. In 
addition inserting the sample in a cavity makes the quan- 
titative analysis more difficult. In Ref ^13i] . propagation is 
investigated in Pr'^+:YSO. The laser is properly stabilized 
but an unexpected low transmission factor is observed. 

In the present paper our aim is twofold. First we pro- 
pose an original solution of the MB equation that analyt- 
ically accounts for the instantaneous part of the atomic 
response. This fast component spreads over a broad spec- 
tral interval and complicates the numerical solution when 
the absorption line is strongly inhomogeneously broad- 
ened. Second, we present a detailed experimental inves- 
tigation of coherent pulse propagation in a large optical 
density Tm'^"*':YAG crystal at liquid helium temperature. 
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Our work is related to topical research on quantum stor- 
age in REIC. The storage protocols generally involve a 
preparation ste p th at can require sophisticated excitation 
procedures [H, 111, [13] ■ We presently concentrate on sin- 
gle rectangular pulse propagation, as a preliminary step 
to the investigation of more complex pulse sequences [T^ . 
In the context of quantum storage, we are more inter- 
ested in the spectral and spatial distribution of atomic 
excitation than in SIT by itself. The 7r-pulse case de- 
serves special attention since all the resonantly excited 
ions should undergo 7r-pulse excitation, wherever they 
sit in the sample depth. 

The paper is organized as follows. In Section 2 we de- 
tail the solution of the MB equation. In section 3 we 
discuss the features of rectangular pulse simulated prop- 
agation. The setup is described in Section 4. We present 
and discuss the experimental data in section 5. 



II. THE INSTANTANEOUS RESPONSE 

In this section we rely on physical arguments in order 
to simplify the computation procedure. Specifically we 
show that one can simply account for the contribution 
from atoms far from resonance. This way the sum over 
the inhomogeneous width can be limited to a narrower 
spectral interval. 

Let us consider a spatially and spectrally uniform dis- 
tribution of two- level atoms. Spectral uniformity means 
that the inhomogeneous width of the optical transition 
is infinite. Let a monochromatic plane wave be directed 
to the sample. At the input side z = the electric field 
E(z,t) reads as: 



E(0,<) = A(0,i) cos(wLi + 



(1) 



where (j) is time- independent. Within the frame of 
the slowly varying amplitude (SVA) and rotating wave 
(RWA) approximations, the Maxwell Bloch equations 
reads as: 



dtu{uJab; Z, t) = -Av(Wafc; Z, t) 

dtv{uJab; Z, t) = -ri(z, t)w{uJab; Z, t) + Au(Wa6; Z, t) 
dt-w{uJab'i Z,t) ^ Q{z,t)v{uJab; Z,t) 

(2) 

where A — Lo^b — and where ri(z,t), ^{LJab', z^t), 
v(ujab]z,t) and w(cjafc;-2,i) respectively represent the 
Rabi frequency and the components of the Bloch vector 
li{ujabT z,t). Let L represent the sample thickness. To 
further simplify the problem, we assume that the incom- 
ing pulse evolution characteristic time, i.e. the inverse 
spectral width 1/Ap, is much longer than L/c. There- 
fore we can neglect the time derivative dt^{z, t)/c in the 
wave equation. 

To solve the equations numerically, one would like to 
get the spatial distribution Cl{z,t) at time t in terms of 



r2(0,t), the boundary value at the sample input, and of 
Ji(ujab', z,t — t), where rAp << 1. Then one would solve 
the Bloch equation at time t, starting with Ji{ujab'-, z, t — r) 
and n{z,t — t). It should be stressed that the material 
response at time t cannot be simply expressed in terms 
of the Bloch vector at t — r. Indeed, due to the infinite 
inhomogeneous bandwidth, the material gives rise to an 
instantaneous response. Let us formally solve the Bloch 
equation in the following way: 



v(wQb; z, t) = v(tJa6; z,t~T) cos(At) 
+VL{uJab; z,t -t) sin(AT) 



ft{z, t — T')w{LLJab', Z,t ~ t') C0s(Ar')(iT' 



(3) 



The first two terms on the right-hand side corresponds 
to the free precession of the Bloch vector from t — t to t. 
The last term reflects the coupling to the field during this 
time interval. Initially w(wa6;^,^o) = ~1, since all the 
atoms sit in the ground state. As the pulse propagates 
through the sample, w{ujab', z,t) + l departs from over a 
spectral interval given by the pulse width Ap. However, 
according to Eq. [31 the driving field uniformly excites the 
atoms all over the inhomogeneous width since the weight 
factor w(ujab] z,to) is initially close to —1 everywhere. 
Summing over an infinite inhomogeneous width clearly 
raises a numerical computation issue. An infinite width 
contribution also reflects an instantaneous response fea- 
ture that we intend to express analytically, thus relaxing 
the computation problem. Substituting v{ujab',z,t) into 
the wave equation one obtains: 

d^niz, ^) = J di^ab ^{z, t - t') cos(AT')dT' 

+ ^ / ^{z,t^T')[w[L0ab\Z,t) + l]cOs(AT')dT' 



~2tT ' '^^'^byfreei^^ab] Z,t - T,t) 



where: 

Vfreeil^ab;Z,t~T,t) = 

v{ujab; z,t~T) cos(At) -I- u(cja{,; z,t-T) sin(AT) 



(4) 
(5) 



The instantaneous response contribution has been iso- 
lated in the first term on the right-hand side. This eas- 
ily reduces to —-^^{z.t). The second term, of order 
— ■|ApTl7(z, t), can be neglected since Apr << 1. Fi- 
nally the wave equation solution reads as: 



Vt{z,t) = r2(0,i)e- 



2^ 7o 



dz'e 2(2 2) I duJab^free{^ab;z' ,t - T,t) 



(6) 



In the small area limit the second term on the right-hand 
side can be neglected. The equation then reduces to the 
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Bouguer law of absorption [l0|. As a starting point 
for numerical computation, the equation conveniently ex- 
presses ri(z, t) in terms of $7(0, t) and of B{u!ab', z,t — t). 
In order to get B{ujab', z,t) in terms of B(ijJaft;z,i — r) 
and ^{z, t — t) one complements Eq. [3] with the integral 
solutions for u(ti'ab; z, t) and w{ujab', z, t): 

u{uJab;Z,t) ^VLfree{i^ab;Z,t-T,t) 

+ / n{z,t - T')w{uJab;z,t - T')sm{AT')dT' 
Jq 

^{^ab]Z,t) =w{uJab;Z,t~T) 



+ / n{z,t - T')viuJab; Z,t - T')dT' 

Jo 

where: 

U/reel^^afc; Z,t-T,t) = 

u{LUab; z,t- t) cos(Ar) - v{ujab; z,t~ t) sin(AT) 



(7) 



(8) 



According to Eq. [71 the second term on the right-hand 
side of Eq. [5] still apparently contains far-from- resonance 
contributions. Indeed \i{ijJab] z,t) and \{ujab]z,t) are 
built from wiujab] z^ t) all over the infinite inhomogeneous 
width. To show the instantaneous response is actually 
contained in the first term, let us define the new vari- 
ables: 



\]{LOab\ Z,t) = 

n{ujab;z,t) + / n{z,t-T')sm{AT')dT' 
Jo 

\{uJab; Z,t) = 

/•OO 

Y{LOab]Z,t)- / r2(z,i-T')c0s(AT')dT' 
^0 



(9) 



Those new variables vanish far from resonance with the 
driving field. Then one easily verifies that Eq. [5] is left 
unchanged if n{LOab] z,t — t) and vicuab', z,t — t) are re- 
spectively replaced by \J{uJab', z,t — T) and Y{ujab', z,t — T). 



III. PULSE DISTORSION AND AREA 
THEOREM 

According to the McCall and Hahn theorem the 
transmitted pulse area, Aout, can be expressed in terms 
of the incoming pulse, Ai^, as: 



tan(Aout/2) 



tan(Ai„/2) 



(10) 



where the pulse area is defined a,s A — J il{t)dt. The cor- 
responding variation of the transmission factor Aout/^in 
as a function of Ai^ is displayed in Fig. [T] The shape 
distorsion of an incoming rectangular pulse is presented 
in Fig. [5] for different pulse area values. The box la- 
bels help to locate the propagation regime in FiglTJ The 
distorsion can be understood in the light of the area the- 
orem and of energy conservation. Four different regions 
can be identified in Fig[TJ In region (I) nearly no light is 
emitted after the incoming pulse extinction. The pulse 
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FIG. 1: area theorem: pulse area transmission factor Aout /Ain 
as a function of the incoming pulse area Ain. Labels (a) to 
(/) refer to the different input-area conditions considered in 
Fig.H 



simply obeys the Bouguer law and is hardly distorted. 
In region (II) , located around tt area, the pulse stretches 
in order to comply with two contradictory prescriptions. 
On the one hand the energy absorption increases, since 
all the resonantly excited atoms are promoted to the up- 
per level, at any depth within the sample. On the other 
hand the area transmission factor grows larger that unity. 
The region (III) is centered on 2tt area. The transmission 
factor is still close to unity but less energy is absorbed 
since the resonantly excited atoms are returned to their 
ground state. In this region the outgoing pulse duration 
reduces, getting closer to that of the incoming pulse. In 
this region the pulse tends to the expected soliton hy- 
perbolic secant shape. Finally the pulse area undergoes 
another increase in region (IV), around Stt area. As in 
region (II), this area increase confiicts with the energy 
depletion by the resonantly excited atoms that are left 
in the upper level. This apparent contradiction is solved 
by the emergence of a stretched secondary component 
in the tail of the main pulse. However, the Stt pulse is 
less distorted than the tt pulse by the increase of energy 
transfer to the atoms. Indeed, to excite the same number 
of atoms, a 37r pulse consumes a 9 times smaller fraction 
of the available energy than a tt pulse. 

The rapid growing of the long tails, in both the tt and 
Stt regions, provides us with key features for experimental 
data analysis. 



IV. EXPERIMENTAL SETUP 

The 0.5% at. Tm3+:YAG crystal is cooled down to 
= 2.2K in a helium bath cryostat. The experiment is 
performed on the ■^114— '^Hg transition at 793nm. The 
sample length is L = 5mm. At the operating tempera- 
ture the opacity is measured to be aL = 5. The crystal 
sides are cut perpendicular to direction [1, -1, 0], and the 
light beam is polarized along direction [1, 1, 1]. This way. 



4 



(a) 


(b) 




0.571 



(c) 

A; = 0.9JI 



1 2 3 4 5 



T3 

3 





(d) 

= 71 


w 

A„= 1.171 








= 1.6 71 


-A 

A A., = 271 




w 

= 371 


-A 

11 A„ = 3.171 


normalized time 





FIG. 2: temporal distorsion of an incoming rectangular pulse 
after propagation through a aL — 5 sample, for different ini- 
tial area values. Time and outgoing pulse amplitude are re- 
spectively normalized to the duration and the amplitude of 
the incoming pulse. 
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FIG. 3: experimental setup. 



three sites out of six interact with the driving field, and 
they do so with the same Rabi frequency. 

The sample is excited by a monochromatic semicon- 
ductor laser (Fig[3]). We carefully control the spectro- 
temporal and spatial properties of the source. The laser 
linewidth is first reduced to less than IkHz by Pound- 
Drever-Hall locking to a high-finesse Fabry-Perot cavity. 
After boost up through a tapered semiconductor ampli- 
fier, the laser beam is precisely temporally shaped by an 
acousto-optic modulator (AOM). The AOM is directly 
driven by a high frequency arbitrary wave form generator 
(Tektronix AWG520), operating at IGs/s, which guaran- 
tees precise amplitude and phase control. The beam is 
then spatially filtered by a 2m-long single mode optical 
fiber. 

The experimental conditions must be consistent with 
the plane wave assumption of the theory. First we make 



the laser beam depth of field (DF) » L. The DF, de- 
fined as twice the Rayleigh range, is given by 27rnwg/A, 
where n — 1.82 stands for the YAG index of refraction 
and where the waist has been adjusted to wq = 25/im. 
Therefore DF= 9mm, which is significantly larger than 
L. In order to put the beam waist at the crystal cen- 
ter, we position the focussing lens at maximum Rabi fre- 
quency, as probed by an optical nutation signal. Then 
we use a pin-hole to select the signal emerging from the 
center of the illuminated spot. The sample is imaged on 
the 50/im-pinhole through a telescope with a x4 mag- 
nification factor. This corresponds to 8% peak-to-peak 
excitation intensity variation over the detected area. The 
signal is detected on an avalanche photodiode. 

The input intensity is monitored by photodetector 
PDcontioi(Figl31)- We again resort to optical nutation 
measurements for calibrating the detector in Rabi fre- 
quency units. Rabi frequency can be precisely deduced 
from the optical nutation signal in the low absorption 
limit only. Because of the large sample opacity at res- 
onance, we have to detune the laser about O.lnm from 
resonance to perform the calibration. The transmitted 
intensity is detected on PDcrystai- We calibrate this de- 
tector with the raising edge of the rectangular pulse. In- 
deed, immediately after switch-on, the transmitted in- 
tensity simply reads as /oexp(— aL), where Iq represents 
the input pulse intensity that is measured on PDcontroi- 



V. RESULTS AND DISCUSSION 

We study the propagation of rectangular-shaped 
pulses. The pulse duration is adjusted to 7/zs, which is 
significantly smaller than the atomic coherence lifetime 
of = 50/is. 

The input area ranges from O.Stt to SAtt. The observed 
temporal profiles are displayed in Fig. [H As expected, 
the pulse is strongly stretched around input tt area, since 
the area has to be conserved despite of large resonant ab- 
sorption. Then, as the input area is increased, the pulse 
shrinks back, which is consistent with reduced absorp- 
tion. Finally, around Stt, the pulse spreads again as reso- 
nant absorption increases, growing a side lobe to simulta- 
neously satisfy area conservation and energy drop. The 
expected output area evolution is also observed. How- 
ever, agreement with theory is only qualitative. 

As seen in Fig. [5] the ratio of the output area to the 
input area lags behind the predicted value, especially in 
the region of tt and Stt incoming pulses. The theoretical 
profile in Fig. [5] accounts for the finite dipole lifetime 
that has been set to 50/is. As compared with Fig. [1] the 
finite coherence lifetime slightly reduces the computed 
transmission factor, especially in the tt and Stt regions 
where the propagating pulse is strongly stretched. How- 
ever this reduction is not sufficient to fit the experiment. 
Finally the excess of absorption reminds us of previously 
reported data p^ . although agreement with theory has 
been improved in our experiment. A possible fault of the 
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A,„ = 0.8ji 




setup is the absence of antireflection coating on the crys- 
tal. As a consequence a reflected field counterpropagates 
through the excited medium, with an amplitude reflec- 
tion coefficient as large as 30%, given the high index of 
refraction of YAG. One may also fear that the plane wave 
assumption breaks down in the 7r-pulse region. Indeed, 
the portions of the Gaussian beam profile correspond- 
ing to areas smaller than tt are absorbed more efficiently 
than the central part of the beam [l2| . Such a stripping 
of the outer portions of the beam tends to reduce the 
waist and to move it to the output of the sample. The 
effective depth of field gets rapidly smaller than the sam- 
ple length, since it varies as the square of the waist. Only 
a 3D Maxwell equation would correctly account for the 
resulting focussing and diffraction effects. 
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FIG. 4: Experimental profile of the transmitted pulse ampli- 
tude (solid line). Input area (Ai„) ranges from O.Svr to SAn. VI. CONCLUSION 
Input pulses (dotted line) are rectangular with 7/is-duration. 
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FIG. 5: experimental and theoretical pulse area transmission 
factor AovLt/Ai-n as a function of the incoming pulse area ^in- 
Input pulses are rectangular with 7/.is-duration. Theory ac- 
counts for a finite coherence fifetime of 50/is 



The propagation of rectangular pulses has been investi- 
gated in monochromatic plane wave conditions. In order 
to consistently account for the very broad inhomogeneous 
bandwidth of the absorption line, the standard MB the- 
ory has been adjusted so that the instantaneous radiative 
response has been described analytically. Experimental 
data qualitatively agree with theory. However a signifi- 
cant quantitative discrepancy subsists, especially in the 
most interesting 7r-pulse region, where the resonantly ex- 
cited Bloch vector should undergo a tt rotation at any 
depth inside the sample [l^]. In the context of quan- 
tum optical storage, the need for efficient preparation of 
larg e optical density sample has been stressed recently 
[2ll | and should stimulate further coherent propagation 
investigation. Extension to three-level A-systems might 
be appropriate. 
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